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Abstract 

We calculate the anomalous dimensions of operators with large global charge J 
in certain strongly coupled conformal field theories in three dimensions, such 
as the 0(2 ) model and the supersymmetric fixed point with a single chiral 
superlield and a W = <E> 3 superpotential. Working in a l/J expansion, we find 
that the large-/ sector of both examples is controlled by a conformally invariant 
effective Lagrangian for a Goldstone boson of the global symmetry. For both 
these theories, we find that the lowest state with charge / is always a scalar 
operator whose dimension A j satisfies the sum rule 

J 2 Aj - (4 + i + A/_i - - i + £) A/ +1 = 0.04067 , 

up to corrections that vanish at large J. The spectrum of low-lying excited 
states is also calculable explcitly: For example, the second-lowest primary 
operator has spin two and dimension Ay + y/3. In the supersymmetric case, the 
dimensions of all half-integer-spin operators lie above the dimensions of the 
integer-spin operators by a gap of order / A. The propagation speeds of the 
Goldstone waves and heavy fermions are and ± \ times the speed of light, 
respectively. These values, including the negative one, arc necessary for the 
consistent realization of the superconformal symmetry at large J. 
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1 Introduction 


Motivation 

Most conformal field theories (CFTs) lack nice limits where they become simple and 
solvable. Some of them are just in the middle of coupling-constant space, and all 
anomalous dimensions are of order 1, and the OPE is just intrinsically complicated, 
and no parameter of the theory can be dialed to a simplifying limit. 

Even in such cases there may sometimes be sectors of the theory where anoma¬ 
lous dimension and OPE coefficients simplify. Many examples of this type are 
known, where the simplifying limit involves taking large quantum numbers, J, under 
rotational or internal global symmetries. 

The simplifications of such limits at large J have played a role in the correspon¬ 
dence between strings in holographic spacetime and single-trace operators in planar 
gauge theory 0 - The large-i limit has also simplified the analysis of the high-spin 
spectrum in CFT [2,3] and the worldsheet theory in the confining regime QQ 
It would be desirable to develop a more general understanding of the simplifica¬ 
tion of CFT at large global symmetry quantum numbers. First of all, in theories that 
are strongly coupled, we should exploit any analytic tools available to gain informa¬ 
tion about the operator spectrum. Second of all, such simplifications in the spectrum 
as have been observed in |3j[7][8] appear to be fundamental to the structure of CFT it¬ 
self as understood through the conformal bootstrap. These simplifications have been 
derived using the abstract rules of the bootstrap rather than any sort of Lagrangian 
formulation^ It is therefore quite intriguing that the form of the asymptotic spec¬ 
trum at large spin is quite reminiscent of the results we derive in this paper for large 
global internal symmetry quantum numbers]^] 


Summary of methods 


In this paper we will illustrate the simplification of CFT at large global symmetry 
quantum numbers with two simple examples of strongly coupled fixed points in 


three dimensions: The critical point of the 0(2) model [ 14[ (or XY model), and the 


jV — 2 superconformal fixed point of the Wess-Zumino model with a single chiral 
superfield and a W — <f> 1 2 3 4 supcrpotcntialQ We treat these theories by quantizing 
them on a 2-sphere of radius R and calculating their operator dimensions via radial 
quantization. 

In the limit of large J, there is a large hierarchy between the radius R of the sphere 
and the length scale set by the charge density, 2 ~ R. We can then consider 
the CFT as a Wilsonian effective action at a renormalization group fixed-point with 


1 See also |5, 61 for earlier attempts at an analysis in the same limit. 

2 For progress on the simple 3D conformal theories using modern bootstrap methods, see for exam- 
pie |9)(I3). 

3 We thank Kallol Sen and Aninda Sinha for discussions and correspondence on this point. 

4 This theory can be shown to flow to a nontrivial fixed point by various means, such as the computation 
of the two-point function of the /^-current in the infrared |T3|[T6| . 
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cutoff A and exploit the cutoff-independence of the dynamics to take 


i«A «VP=Si- < U) 

which implies 7 > 1. In the regime A <C yfp , it turns out that the effective action 
is weakly coupled and under perturbative control, with an expansion in powers of 
p _1 , with numerators given by derivatives of p and powers of the cutoff A. When we 
quantize the theory, the leading approximation to any quantity is given by the leading 
large-p term in the action. Sub-leading corrections are generated by quantum loops 
and by explicit vertices with more negative powers of p in the operators. 

In the Wilsonian action, only non-negative powers of A appear, and these are 
small in the limit ( |1.1| ). Furthermore, the underlying conformal invariance of the 
theory means that the A-dependence cancels in all physical observables, order by 
order in the j expansion. The A-dependent terms are scheme-dependent. They play 
little role in the dynamics of the large-7 theory other than to restore quantum scale 
invariance in the effective Lagrangian description and to cancel the A-dependence 
in amplitudes. We will explain how the A-dependent terms may be calculated algo¬ 
rithmically order by order from the A-independent terms for any given form of the 
cutoff. 


Summary of results 

In the three-dimensional examples we consider, we find that leading large-7 behavior 
of the dimension Aj of the lowest operator with global charge 7 goes as 

Ajocj+i (1.2) 

at large 7. This is true both in the 0(2) model and also in the W = <f> 3 model. 
Note that the dimension in the latter case does not go as A = 7 + 0(J°), as might 
naively have been expected based on supersymmetric considerations. Despite the 
presence of a BPS bound and multiplet-shortening condition in this -sV = 2 theory, 
the operators (j) J do not saturate it, even approximately: The lowest state in the 
large-7 sector is parametrically far above the supersymmetric bound. 

We can also compute the sub-leading terms in the expansion. In both the 0(2) 
model and the W = <f> 3 model, the lowest dimension in the large-7 sector has the 
expansion 

A y = C3 7 3/2 To7 1/2 -0.108451 + 0(7~?) . (1.3) 

2 2 

at large 7. The coefficients c 3 and ci are related to the coefficients of the leading 

2 2 

and first sub-leading terms in the large-7 effective Lagrangian and we do not at 
present know how to calculate these coefficients from first principles. They may 
differ between the supersymmetric and bosonic models. The order 7° term, on the 
other hand, is calculable and common to the two models, and to any other model 
described by the same large-7 universality class. 
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Dimensional analysis and large-J scaling 

We note here that the leading large-7 scaling can be deduced immediately on dimen¬ 
sional grounds without the need for the methods developed in this paper. At large 
charge, the charge density p and energy density Jif are homogeneous and semiclas- 
sical on distance scales between R and p~ 2 . In this range of scales, the two densities 
must obey a local and scale-invariant relationship, which in three dimensions can 
only be of the form 

jT-p+i (1.4) 


at large J, leading to the scaling (1.2). 


The only exceptions to the rule (1.4) in three dimensional CFT are theories that 
have a vacuum manifold of exactly flat directions, such as a free complex scalar or a 
supersymmetric theory with a quantum mechanically supersymmetrically protected 
moduli space. For such theories, the size of the sphere is never irrelevant, because 
in the absence of the conformal coupling term of the scalar fields to the Ricci cur¬ 
vature, the spectrum of the Hamiltonian would collapse and become continuous. 
Such theories do indeed have continuous spectra on flat spatial slices such as T 2 , 
and therefore the Ricci curvature Ric 3 must enter into the leading term in any local 
relationship between the charge density and energy density, and such theories obey 
3%* ~ (Ric 3) + 2 p for the ground state at large J. 

CFT with good spectra on flat slices, on the other hand - by which we mean, 
theories with discrete energy levels when quantized on T 2 - are of a more generic 
character, and include all known interacting non-supersymmetric theories, and even 
many supersymmetric ones. The critical 0(2) model, the W — <t> 3 model, three- 
dimensional super-QED and the CP(tV) model [17], non-supersymmetric Chem- 
Simons theories with matter [TT8l IT9]], and even theories of free fermions all fall 


into this category. The scaling of the operator dimension with global charge in such 

3 

theories is always J 1 -. This scaling is directly visible, for instance, in the theory of a 
free relativistic fermion, where the 7 + 2 is just the ground state energy of a Fermi sur¬ 
face on the sphere. The same scaling also appears in interacting three-dimensional 
theories, such as Chem-Simons theory with matter, at large magnetic flux number]^] 
Despite the automaticity of the leading-order 7 + 2 scaling itself, the power laws 
appearing in the sub-leading large-7 corrections do not follow directly from dimen¬ 
sional analysis, nor does the computability of the coefficients of those corrections 
in perturbation theory. Rather, the structure of the sub-leading corrections follows 
from a renormalisation-group analysis that may depend on the details of the theory. 


Outline 


The plan of the paper is as follows. In Section [2] we discuss the RG flow and 
large-7 perturbativity based on a toy model. In Section [3] we apply the lessons 


We thank Ethan Dyer and Mark Mezei for discussions and correspondence on their results for this 


example and related work by them 17 20 211. 
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we have learned to the full 0(2) model in three dimensions and discuss its large-7 
analysis and study its reduction to Goldstones after integrating out the a field. We 
go on to discuss the classification of operators in the conformal Goldstone action. 
In Section [ 4 ] we move on to the supersymmetric W = <f> 3 model. We perform its 
RG analysis at large 7, decouple the fermions and compute the large-7 expansion of 
operator dimensions using radial quantization. In Section [5j we work out the energy 
of the excited states and see which states are primaries and which are descendants. 
In Section [6| we go on to discuss other models in three and four dimensions. In 
Section[7J we present our conclusions and point out interesting future directions. 

2 RG flow with a dimensionful vev 

We would like to investigate to what extent the operator spectrum simplifies and be¬ 
comes calculable at large 7 in generic CFT, particularly in models for which there are 
no other tools or tunable parameters of the Hamiltonian. A paradigmatic example of 
such a theory is the critical point of the 0(2) model in three dimensions, whose ul- 

traviolet definition is simply a complex scalar field with potential term V = \ (j> | 4 . 

On the one hand one expects the spectrum of this model to be under perturbative 
control for large 7: The loop counting parameter for a process with characteristic 
scale R and initial and final state |7) is y — g 2 R/J. On the other hand, the energy of 
the 7— state on a sphere of radius R scales as J/R and becomes so large as to exceed 
the scale where the ultraviolet physics decouples if we require y -C 1 ,i.e. 7/k> g 2 . 
This would suggest that large-7 perturbativity is not useful for the computation of 
operator dimensions in the interacting CFT. Such a conclusion is too pessimistic. 
We shall see how to compute certain quantities in the large-7 sector in a controlled 
fashion. 

We first turn to the analysis of the renormalization group equation for this sys¬ 
tem, in states of large 7. We will begin with a toy model of the RG flow, that can be 
thought of as the 0(2) model where only the flow of one particular operator is re¬ 
tained. After understanding the general behavior of the toy model, we will return to 
the full 0(2) model and analyze its behavior at large charge density. Then, we will 
solve the full RG equations at the fixed point, in the large-7 expansion. The solution 
in this limit reduces to an effective Lagrangian that is classically scale invariant, plus 
small quantum corrections suppressed by positive powers of the ratio of the cutoff 
to the square root of the charge density. 

This Lagrangian explicitly realizes the conformal symmetry of the underlying 
CFT, even while strongly spontaneously broken by the charge density itself. We then 
quantize this theory on the unit S 2 spatial slice, and calculate the energy. This gives 
the value of the dimensions of operators at large 7 via radial quantization according 
to the state-operator correspondence. The energy of the lowest state at large 7 can be 
calculated in an asymptotic expansion in inverse powers of 7 in terms of our effective 
Lagrangian. Certain terms in the expansion depend on the coefficients in the large-7 
Lagrangian, while other terms are universal and theory-independent altogether. 
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2.1 Toy model of the RG flow 


To understand the effect of large J on the renormalisation group flow, let us write a 
toy model for the renormalization group equation for the effective action of the 0(2) 
model. In this model, we discard renormalizations of operators with derivatives and 
also renormalizations of the mass term and other terms in the potential itself, other 
than the quartic coupling. 


We define this model in the spirit of the toy model in [ 221: We take the structure 
of a true exact Wilsonian RG flow for a real scalar field, truncated to the flow for 
the quartic coupling. We discard the cutoff dependence of the mass term, the ki netic 
term, and all higher-derivative terms. Though this is an uncontrolled approximation 


of the full RG flow, as in the toy model of [22], it illustrates a key behavior of renor¬ 


malisation group flow in a situation where the infrared quantum scale invariance is 
strongly spontaneously broken by the expectation value of a field. 

The ultraviolet Lagrangian will be that of a three-dimensional real scalar with a 
quartic potential: 


2 %, 


1 4 

"2 <!,a > “l2 fl 


( 2 . 1 ) 


Conformal behavior for A Cg 2 


We integrate out modes and lower the cutoff, which is recorded by a A-dependent 
evolution of the effective Lagrangian , of which we keep track of only the renor¬ 
malized quartic coupling g^ eff j(A). At A ~ g 2 , the quantum corrections to the cou¬ 
pling (A) are of the same order as g 2 itself. At that point, the coupling g 2 c|r (A) 

quickly reaches its attractive)^] fixed-point value in units of A, after which point it 
obeys 


ffjeff] ( Al ) _ g[eff]( A2 > 
Al A2 


Al,2 < g 2 , 


( 2 . 2 ) 


or equivalently 

( A ) = hA , A <g 2 , 


(2.3) 


where h is a dimensionless coupling whose numerical value is determined by the 
fixed point equation. 


VE V for the a field as an infrared cutoff 

We can calculate the effective Wilsonian action for any value we like of the modes 
below the cutoff. In particular, we would like to give a nonzero value to the constant 


6 In the full field theory, the fixed point is not fully attractive; one must fine-tune the mass term for a 
in order to stay on the fixed point. In our model, however, we have simply truncated out the quantum 
running of the infrared-unstable coupling m 2 a 2 along with all the higher-derivative couplings we 
don’t want to keep track of. 
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mode of the scalar field a. We assume this value is far below the scale set by the 
ultraviolet scale, so a <^g. This expectation value induces a mass 

m l =4ff] a2 (24) 

for the ^-fluctuations. When we integrate down to momentum modes comparable 
to or below m a , the evolution of the couplings stops, because modes of a no longer 
make large contributions to the running when they are below the gap. Fluctuations 
of such modes renormalize the couplings only with a suppression by positive powers 
of —. Therefore, at approximately 


Amass — 


= m, 


a 


(2.5) 


the evolution of g*- eff ^ comes to a halt, 
infrared value of the coupling satisfies 


Combining equations (2.4) and (2.5), the 


k[eff] (Amass) — h~ a" , (2.6) 

Amass ~ m a = hoT . (2.7) 


Then the effective potential below A mass , which is V[ e ff] ( a ) = j 28 [eS\ (A ma ss) « 4 be¬ 
comes 

V m (a) = ■ ( 2 - 8 ) 

This form for the low-energy effective potential could have been inferred just on 
the basis of dimensional analysis, but it is instructive to see how the sextic potential 
emerges from the structure of the renormalization group. 

For A < A mass = her. the effective action asymptotes to a constant as A —y 0. De¬ 
spite this, the theory is still conformally invariant; the conformal invariance is simply 
broken spontaneously by the expectation value of a. Since quantum corrections to 
the RG flow are small in this regime, the effective action should be classically scale 
invariant to first approximation, with small quantum corrections. In quantum terms, 
classical scale invariance is broken by operators appearing with coefficients propor¬ 
tional to positive powers the cutoff. This effective action has a smooth limit as A is 
taken to zero, since there are no long-range degrees of freedom to generate infrared 
singularities. 


2.2 Coupling to another sector 

Dynamical source for the a-field 

Now suppose further that the vev of a is not fixed arbitrarily, but set by a source 
term B(x)a 2 , in which case the vacuum expectation value of a goes as R + ?, and for 
A <C A mass ~ B + 2, the effective action is classically scale-invariant as a functional 
of B and its derivatives, with B assigned a scaling dimension of 2, and quantum 
corrections suppressed by positive powers of -j=. 
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We can go further and let B(x) itself be a dynamical object, constructed as an 
operator in a local effective field theory coupled to the a sector. This theory must 
also be regularized and renormalized in its own right. However if the sector in which 
the B operator lives also has trivial infrared dynamics - either because it is gapped or 
because it becomes free and massless at low energies - then once again, the effective 
Lagrangian will have only positive powers of A at the lowest energies. 

Integrating out the afield 

For purposes of examining physics of the B sector at scales below A mas s, we can 
go further still and integrate out the a field altogether. We are then left with some 
effective theory of the B sector. Here, powers of B can appear in the denominators 
of terms in in the effective theory, because B sets the mass scale for the a field. 
However the derivatives of B only appear as polynomials in the numerator, in the 
expansion where B is varying slowly compared to the length scale \B\ ~2 itself. 

If the coupled dynamics of a and B are conformal below the UV scale g 2 , then 
the effective action in the /Lscclor must itself be invariant under the conformal sym¬ 
metry at scales below the mass of a. This theory is under control when the conformal 
symmetry is broken spontaneously, since s/B will always be comparable to the mass 

of a. The effective action will contain terms such as |. 6 | + 2 , ^ ^ , et cetera. 

1*1 3 1 * 1 2 

Quantum mechanics of the effective B-theory 

We can then quantize this effective theory in an expansion where the rate of variation 
of B is small compared to the cutoff A, which in turn is small compared to the scale 
s/B ~ A mass ~ a~\ 

^<A Cs/B. (2.9) 

In this expansion, quantum corrections to amplitudes are suppressed by p/VB 
and A/ \[B. As a result, the effective theory of the B sector is conformally invariant, 
not only quantum mechanically, but approximately classically as well. Corrections 
to classical scale invariance in the effective Lagrangian are given by terms with pos¬ 
itive powers of A/ \[B. In the case where the full theory is conformally invariant, the 
quantum mechanical, cutoff-dependent terms can be determined systematically in 
terms of the classical, cutoff-independent terms through the renormalization group 
equation. 

Let the effective Lagrangian for the B sector be a sum of terms 

^ {A] [B}=^ + ^ , ( 2 . 10 ) 

where the terms are classically scale-invariant, and the terms 2 ^j u 1 are not. 
That isQ 

^A, U) = E A 3 - t 4 q ” ) . (2.11) 

A f ^3 

7 In fact, we can allow a slightly more general form of the expansion, with logarithms ln(Z?/A 2 ) on 
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where all terms in 2 z?have classical scaling dimension A c equal to 3, and „^ qu) 
has classical scaling dimension equal to A c . 

If 2z? is the renormalization of the Lagrangian as we integrate out a shell 
of modes of infinitesimal thickness 8 A, then the equation for the theory to live at a 
fixed point is simply 


5 ( rg ) 

A ^= 


E (Ac-3) A 3 " 

Ac 7^3 


■ c/>{ qu) 


( 2 . 12 ) 


Now, we are working in the limit (|2.9|), so can expand both sides in powers of /=. 


By our assumption, the Lagrangiah tor the B sector is trivial in the infrared (either 
gapped or else free and massless), so the left-hand side has only terms with A c < 3. 

This does not mean, however, that the right-hand side contains only a finite num¬ 
ber of terms. Since we are working in an effective theory where the operators are 
allowed to contain powers of B in the denominator, the classical scaling dimensions 
of the operators appearing in the A-dependent terms can be negative: The set of 
allowed values of A c is bounded above by +3, but unbounded below. 

For instance, the RG flow of the classical piece of the Lagrangian may take the 
form 

<5( RG ) „ . , „ A 5 ( dB) 2 


A—— ^K l A 3 + K 2 — + K 3 A 5 
8 A B 

Then the leading pieces of the quantum action must be 


B 4 


■ + ■ 


(2.13) 


<v>(q u ) _ 

-HA] “ 


— K\ A 3 --K 2 — 
3 5 B 


1 

5 


K 2 A 5 



(2.14) 


Thus we can "bootstrap" the coefficients of the operators in the A-dependent 
quantum terms 2zf( qu ) algorithmically from the form of the classically scale-invariant 
Lagrangian (cl ) . though 2z?( cl ) may itself have unknown coefficients. 


2.3 Lessons for large-/ perturbativity 

Though the discussion at the level of generality in this section may be unfamiliar, 
many examples of this type of dynamics are already well known. The most familiar 
examples occur in supersymmetric theories with moduli spaces of exactly super- 
symmetric vacua. In these theories, the role of the B sector is played by the massless 
moduli fields and unbroken gauge fields themselves, while the role of the a sector is 
played by the massive, spontaneously broken gauge fields and the scalars that gain 
a potential through F- and D-tcrms. 

The toy model above contains the essence of the large-7 analysis of the critical 
0(2) model in three dimensions and that of its supersymmetric cousin, the super- 
conformal IT = T> 3 model. In these theories, the role of the a field is played by the 

top of powers of A. The logarithms are physically quite interesting, but somewhat beside the point of 
the discussion here, so we do not emphasize them. The only potentially important logarithm would 
be a term in the expansion with cutoff dependence A 0 (ln(B/A 2 )) k for k > 0. But we have assumed 
the theory of the B sector to be trivial as A —> 0, which does not allow for such terms. 
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magnitude of the complex scalar (j), and the B operator is constructed from the gra¬ 
dient squared of the phase variable %■> that is, B(x) = \d%\ 2 , which is proportional to 
the charge density of the system. 

Our large-7 descriptions of the conformal 0(2) model and W = <f> 3 model will 
be effective theories for the phase variable x after integrating out the magnitude 
a = \(f>\ of the complex scalar. In the supersymmetric W — <f> 3 theory, the fermions 
will also acquire a gap of order m a ~ y/f = \d%\, and be integrated out, leaving us 
with the effective theory of x alone. 

Keeping in mind the toy example of this section will guide us through the theory- 
dependent details of the renormalisation-group analysis of the critical 0(2) model 
and superconformal W = <t> 3 model. 

3 Bosonic 0 ( 2 ) model in 3D 

3.1 0(2) model and its large-/ analysis 

Now we apply the lessons of the toy model to the full 0(2) model. The UV La- 
grangian of the system, which is stable and renormalizable, is of the form 

^uv = -^d M 0-c|0| 2 -g 2 |4>| 4 , (3.1) 

where c is fine-tuned in order for the Lagrangian to flow to a conformal IR fixed 
point [14]: 

UV theory RG flow > IR conformal fixed point. (3.2) 

In order for us to make conceptual contact with the toy problem of Section[2] we 
parametrize (j) as 

(j) — ae' x , (3.3) 

where a G M + and x is 2^-periodic. 

The full exact Wilsonian RG equations are complicated. They simplify quite a 
bit in the large-7 limit, but not until the cutoff is lowered to ■ s /p. The charge density 
p cannot affect the flow between the scale Auv = g 2 and the scale A = x fp. Between 
these scales, the renormalisation group flow is unaffected by the vev. In this range, 
the equations contain all the complexity of the full Landau-Ginzburg theory in its 
strongly coupled regime, including its flow to the conformal fixed point. This range 
of the RG flow can be treated only numerically in a standard treatment, and is not 
simplified by large 7 at all. It is only when we reach the scale A = v /p that the RG 
equations simplify. 

We find, however, that if we use as input the fact that this theory flows to a 
fixed point, we can strongly constrain the form of its Wilsonian effective action in 
the regime of large charge density. When p 1 2 is much larger than the cutoff, the 
magnitude field a decouples as in the toy model of the last section. Let us now study 
the effective action for the complex scalar in the regime 

ACa 2 < g 2 . (3.4) 


10 



Approximate Classical Scale Invariance 


In this limit we write the Lagrangian asj^] 

J?ir = — ^(d^a ) 2 — f(a)(d ll x ) 2 — V(a) + (higher derivative terms). (3.5) 

At the IR fixed point, the Lagrangian is approximately classically scale-invariant, 
with corrections to classical scale invariance that go as A/a 2 . The leading La¬ 
grangian density must have mass dimension 3, and the fields have dimensions 

a [mass] 1 / 2 , % [mass] 0 . (3.6) 

The former is fixed by the normalization convention for the kinetic term of a, and 
the latter is fixed by the dimensionless periodic identification of %. 

It follows that the functions in the Lagrangian have to scale as 

V(a)<*a 6 , f(a)°za 2 . (3.7) 

Therefore the Lagrangian at the IR fixed point can be written as 


tfR = A(d lt a) 2 ~ 1 -Ka 1 (3 llX ) 2 



+ (Ricci coupling) 


+ (higher derivative terms), 


(3.8) 


where k and h 2 are numerical constants. 

Note that the equilibrium value of a lies at 

‘• i = -%(*X) 2 = + 2 1 $\3x\ 2 - (3.9) 

(Here, we have used the fact that the gradient of % is timelike around a state of large 
charge density.) The charge density p = Ka 2 % is 


P = 


Kh 2 


and its integral on the sphere is 


J = 


f d 2 xp = 4nR" 

Js 2 


Kh 2 


In terms of the charge density p, our double hierarchy (3.4) becomes 

A < 


(3.10) 


(3.11) 

(3.12) 


Since we are working on a sphere of radius R, in order for the effective action to 
be consistent, A must be parametrically larger than the IR cutoff M m — l/R. This 
means that A is in the regime 

1 R <<A<< S^ <<gl <3 ' 13) 

which can only be consistent if the charge is large: 


7>1. (3.14) 

The normalization condition for the field has been chosen so that the kinetic term of a is unit and 
canonical. 
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Renormalization Group Analysis 


Let us now analyze the RG equations for the effective action (3.5), using the property 
that the underlying theory is conformally invariant. In particular, let us show that 
the corrections to classical scale invariance are small, and that the cutoff-dependent 
quantum terms can be derived algorithmically from the cutoff-independent classical 
terms by the renormalization group equation. 

Now let us prove that the classical scaling above is correct for large enough 
charge density p. For this purpose, we want to estimate the change of f(a) and 
V(a) in (3.5) with respect to the change in the momentum cut-off A. 

We first Taylor expand around the minimum of a, which will be denoted a^. We 
also let a = a — ao'. 


-2ht = - 

- V (flo) - av'(a 0 ) - ^aV"(a 0 ) H - 

-/(«o)(Wx) 2 -«/(ao)(c>„x ) 2 - UYMtfuX ) 2 + ■■■ 


+ (higher derivative terms) 
Here the Feynman rules are as follows : 

a ■■ 

k 


X - 

k 




k 2 + V"(a 0 ) 

1 

2 f(a 0 )k 2 

(3.16) 

(3.17) 

f\a 0 )(ki-k 2 ) 

(3.18) 

\f"( a o)(k\ -k 2 ) 

(3.19) 


(3.20) 


The double lines represent a-propagators and the single lines represent % propaga¬ 
tors. Apart from the ones shown above, we have one more diagram representing the 
source term aV'(ao): 

(g) = V'(ao) (3.21) 

First we renormalize the one-point function. As we have defined the effective 
potential is at its minimum when a — 0, the full one-point function becomes zero. 
Therefore, we shall not worry about one-point functions anymore. 
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Diagrams that contribute to the renormalization of f(a ) at 1-loop are 

P 


(a) 



(3.22) 



We now calculate these diagrams. We will take the momentum cut-off A and inte¬ 
grate the loop momentum over a shell with radius A and thickness 8 A. 

(a) The contribution from diagram (a), S\f(cio), is given by 


Sif(ao) oc — 


A+8A 



A 


fc 2 /"(ao) 
p 2 + V"{a 0 ) 


f"{a 0 )A 2 8 A f"(ao)A 2 8 A 
A 2 + V"[oq) ~ V"(a 0 ) 


(3.24) 


The last approximation holds when p/A 2 is large enough: ciq oc p 1 / 4 becomes 
large enough to satisfy V"(ao) S> A 2 in our limit p 3> A 2 . 

(b) The contribution from diagram (b), 82 /{a 0 ), is given by 


$»/(« 0 ) « ^ 


A+8A 


d 3 p 


1 


(f'(ao)) 2 (k-p) : 


A 

A+5A 


{k-p ) 2 + V"(a 0 ) f{a 0 )p 2 

cos(0) = l 


I d \p\\p\ 2 I 

A cos(0)=—1 

(/'( ao )) 2 A 2 «A 


J(cos(0)) 


(/ / (q o )) 2 )A-cos 2 (0) (3.25) 

V"{a 0 )f{a 0 )A 2 


V"{a 0 )f{a 0 ) 

Again the approximation holds when p is large enough compared to A. 

Summing over the two contributions above, we find the RG equation below for 
the quantum scaling (assuming that f(a) °c a n )\ 

/(«0) A 3 


4 


/(« 0) 


V"(a 0 ) 


(3.26) 


lq 0 

We can compare this with its classical scaling, which is of course just proportional 
to f(ao) itself. Therefore the ratio of the quantum scaling over the classical one is 


A -A 
Sa 


/(« 0) 


/(« 0) 


(3.27) 
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This is small when p is large enough, i.e. when a 2 ^ A. The analysis for V(ao) 
goes parallel to the analysis above: 


A 


_ 8 _ 

8A 


V(a 0 ) 


V(ao) 



(3.28) 


Just as in the discussion in section [272] we can derive all A-dependent coefficients in 
the effective Lagrangian from the classically scale-invariant ones, through the fixed 
point condition in the renormalization group equation. 

In the limit A —» 0 with the field configuration and in particular a f 0 held fixed, 
the Wilsonian action is finite. This is clear, because the only possible singularities 
as A —» 0 would come from integrals at low momentum. Since a is gapped and X 
is infrared-free in the absence of interactions with a , there are no diagrams that can 
generate such singularities. 


3.2 Reduction to Goldstones 


Classical elimination of a 


Starting with (3.8), we now integrate out the a field and examine the effective action 
for the x field. The mass term of the a field is proportional to \dx\, so in the limit 


\dx\ >A we can set the a field classically to its equilibrium value (3.9), with the 
quantum co rrect ions from its fluctuations being proportional to positive powers of 


Using (|3.8|) we have 


££ = b x |(9^| 3 + (lower order in \dx\) , 


(3.29) 


where] <9^ | = [-{dx ) 2 ] 2 and b x = ^^~. 


Terms of lower order in dx 

The higher-derivative terms come from the kinetic term for a , as well as from higher- 


derivative terms in the effective action (3.8). The kinetic term for a generates higher- 
derivative terms such as . This term comes along with a coupling to the 


\9X\ 

Ricci scalar, proportional to Ric 3 \dx\- The coefficients of these two terms are not 
independent; Weyl invariance relates them to one another. 

There is an infinite series of such terms in the effective action, with arbitrary 
derivatives of x in the numerator, and only \dx\ occurring to negative powers. These 
operators are classically scale-invariant, and are arranged hierarchically in terms of 
the number of powers of dx that occur in them. This is the natural organization 
of operators when we compute observables in a state of large and approximately 
constant density, 

— < A < J~p . (3.30) 

P 
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In Subsection 3.4[ we will classify the first few leading operators in this expansion 
of ££ j^] First, we briefly analyze the quantum properties of the conformal Goldstone 
theory. 


3.3 Conformal Goldstones at the quantum level 


At first sight, the action ( 3.29| ) looks quite strange. The form of the Lagrangian 
density is singular and infinitely strongly coupled as an action expanded around the 
origin % — 0. But the action ( |3.29[ ) is not meant to be used there: It is only a first 
term in an infinite series of terms with higher powers of \d%\ in the denominator. 

Rather, the action is meant to be expanded around a background dxo °c ^/po 
of constant charge density, and the field % quantized in small fluctuations whose 


wavelength is long compared to p 0 2 . In this limit, the Lagrangian is fully under 
perturbative control. 

To see this, break up % into a background and fluctuation, X = Xo + Zfluc> where 
X has constant gradient. The expansion around a fixed background Xo takes the 

formal 


& = b x || 3 + \dXo\ [quadratic in define] + [cubic in define] +6>( ^^ ) , 

(3.31) 

So the quantum mechanical fluctuations of Xfluc are of order \dXo\~ 1 , or equiv¬ 
alently order p _ L It is sometimes convenient to work in terms of a scalar field 
X = |P^o| + 2 Zfiuc with unit-normalized ki netic term. In terms of the hatted Gold- 
stone, we have 


& = b x \ d Xo\ 3 + [quadratic in d%] + \d%o\ 2 [cubic in d%] +^(^r) • ( 3 -32) 

The quantum fluctuations of the canonically normalised field % are of order 1, and 
all the p-dependence comes through coefficients of the vertices which are powers 
of dxo — 0(p + 2). The effective coupling constant that suppresses each successive 
quantum loop is therefore \dXo\ 3 ^ p 2 . 

There is a simple rule to keep track of the quantum-mechanical p-scaling of a 
given term, without explicitly breaking up the x field into background and fluctua¬ 
tion. Since the background Xo has constant gradient, then the term dx is dominated 
by the background piece dxo, while the terms d k x lor k >2 are purely fluctuation 


9 This is in the spirit of the functional RG flow (see j23]| and references therein). Note that our analysis 
is novel because we are not making approximations such as large N or e expansion, but we hold the 
dimensionality of space and the field content constant. 

10 We have dropped the term linear in the fluctuations, which is a total derivative. This term is important 
in that it contributes to the form of the expression for the conjugate momentum to x and therefore 
the charge density, in the small-fluctuation expansion. However it does not affect the equations of 
motion or the computation of corrections, so we omit it in the present discussion. 
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pieces, since d k Xo — 0. So the rules are 


dX ~ dxo = 0(p- 


d k x = d k Xnuc = \dX()\ ~ 2 d k x = 0(p 


(3.33) 

(3.34) 


These rules will allow us to make quick work of the enumeration of the first sub¬ 
leading operators in the Goldstone effective action, in Section [T4j 

Cutoff-dependent terms 

The effective Lagrangian for X is classically scale-invariant in the large-p limit, and 
also fully quantum-mechanically scale invariant by virtue of our input assumption 
that it is the effective theory of an exactly conformally invariant theory. As a re¬ 
sult, one needs in general to include cutoff-dependent terms to restore conformal 
invariance at the quantum mechanical level. 


Just as in the case of the toy model of Section [2T2| we can compute the quantum 
terms algorithmically in a A 3 /|<9^| 3 expansion, from the starting point of the A- 
independent, classically scale-invariant terms in the Lagrangian. The form of these 
can be seen by taking the cutoff-dependent terms in the (j) system and reducing those 
to Goldstone terms by integrating out a — \(p\. Alternately, a renormalization group 
equation for the x action determines the cutoff-dependent terms directly, in terms of 
the canonically scale-invariant terms in the Goldstone action. The quantum pieces 

3 A 3 (d 2 x) 2 


will consist of terms such as A J 


\dx\ 4 


The cutoff-dependent terms are of course scheme-dependent and not very in¬ 
teresting in themselves, however important they may be as a point of principle for 
establishing the self-consistency of our treatment of the theory at large J. For low- 
order calculations, such as the one we will perform in Section |5(2} an analytic renor¬ 
malization scheme such as dimensional regularization or ((-function regularization, 
is far more convenient, as it subtracts the counterterms automatically and we never 
have to calculate the explicit form of the quantum terms in the action. In Section [53j 
we will indeed use ((-function regularization, and therefore we will not pursue the 
calculation of the cutoff-dependent terms further. 


3.4 Classification of operators in the conformal Goldstone action 

Let us now classify all possible scale-invariant operators appearing in the conformal 
Goldstone action, according to their p-scaling. In order to understand the p-scaling 


of a given operator at the quantum level, we use the simple p-counting rules ( 3.33 ), 
( |3.34 ) derived in Section [33] 

We are going to write down scale-invariant operators at order p a -°. The scalings 
to remember are dx ~ p 1//2 and d ... dx ~ p -1 / 4 . We retain only scalar operators of 
scaling dimension 3, including curvatures of the background metric in our counting. 
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3 

Order p + 2 

There is only one operator of order p + 2 : The leading-order Lagrangian 

^3 = \dx \ 3 ■ (3.35) 

This term is conformally invariant as well as scale-invariant, as an operator in the 
Lagrangian density. 

Order p + 2 

At this order, there are two operators and they are both curvature-dependent. The 
first is 


=Ric 3 |<9z| , (3.36) 

where by Ric 3 we mean the three-dimensional Ricci scalar. This term is scale- 
invariant as a term in the action but not conformally invariant. This is related to the 
fact that the Ricci scalar transforms inhomogeneously under a Weyl transformation 
of the background metric g„ —> exp{2u(jt)}g... To restore Weyl-invariance and 
conformal invariance, this term must be completed with an operator of the form 
in the Weyl-invariant combination 


= Ric 3 \dx\ +2 


(d\d x \) 2 

\dx\ 


Note that the p-scaling of the completing term is negative. 
There is also the operator 


RM^X)(d v X) 

\dx\ 


(3.37) 


(3.38) 


This operator is of order p + 2 for a generic background metric. However we will be 
mostly interested in the case of an unwarped product metric of the form (time) x 
(spatial slice), and a spatially homogeneous solution Xo- For such a metric, 
the only non-vanishing components of Ru v are the spatial components, which have 
vanishing contraction with the background gradient d^Xo- Therefore this operator 
contributes only through fluctuation terms. Replacing each of the two contracted 
dxo’s with a dxtiuc = \dX 0\~ 1 d%, the scaling of the operator comes down to p~ J , 
which is below 0 and so we discard it. 


Order p + 4 


We obtain two possible terms at order p + 4 ; 


.t/4 




_ 1/4 

- rs -' o 

\d%\ P 


(3.39) 

(3.40) 
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These terms are not linearly independent as operators; one combination vanishes 
by virtue of the equations of motion. The remaining linear combination is non¬ 
vanishing as an operator, but does not appear in the effective Lagrangian, because it 
is odd under the charge-conjugation symmetry of the theory: 

X~>~X- (3-41) 

Hence there are no terms of order p + ? appearing in the Lagrangian. 

Order p° 

There are no operators at order p°. This is a key fact so let us establish it formally. 
Since first derivatives of % scale as p + 2 and higher-derivatives of % scale as p - ?, the 
number of %’s with higher derivatives on them must be even for an operator to scale 
as a half-integer power of p. Dimension-3 operators with <9%’s only are of the form 
(dx) 3 k ^k, where % is a geometric invariant of scaling dimension k, constructed 
from background curvatures. For k — 3 there are no such operators p] 

The next possibility would be terms with two or more %’s in the numerator with 
higher derivatives acting on them. The total dimension of the numerator would be at 
least four, and have p scaling bounded above by p _ 2 . The dimension of the operator 
could then only be brought down to 3 by putting one or more factors of <9% in the 
denominator, which would drive the total p-scaling of the operator to p 1 or less. 

We conclude there are no operators of order p°. This fact has physical signifi¬ 
cance, and we shall see in Section [ 5 ] that the J° term in the expansion of the operator 
dimension at large J, is calculable and universal as a result. 


4 The Supersymmetric W = O 3 Model 


4.1 The model 

Consider the ,jV = 2 supersymmetric theory in three space-time dimensions, with 
a single chiral superfield, Kahler potential K — and superpotential W — 3<f>\ 
This theory is known to flow to an interacting superconformal fixed point, as can 
be shown, e.g., by the techniques of [ fT5| T6| . The iCcharge and dimension of <£> at 
the fixed point are equal to +|. This theory has no marginal deformations and no 
small parameter of any kind. Nonetheless we would like to analyze its spectrum 
of operator dimensions. It is equipped with a continuous global symmetry, namely 
the /^-symmetry itself, and one can inquire again about the dimension of the lowest 
operator |/) of ^-charge /. (We define the charge here to be | times the P-chargc, so 
that <j> has 0-charge +1 and the supercurrent Q a has <j>-charge — When referring 
to the "charge" without specifying, we will mean the charge normalised as 0-charge 
rather than P-chargc.) 


11 The gravitational Chern-Simons term, which has dimension 3, is not gauge-invariant when multiplied 
by other operators such as dx- 
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4.2 RG flow in the W = <E> 3 model 


Following both our toy model and the example of the 0(2) model in Section [3j we 
expect large-7 perturbativity to mean that the flow of the Lagrangian with log(A) 
is parametrically suppressed relative to the original Lagrangian. This is indeed the 
case; the calculation is mostly parallel to that case and so we suppress most details. 

In our normalization convention for the fields, W is taken to be unit-normalised 
(and thus has scaling [mass] 2 ). This is different from the standard convention in 
which the kinetic term is unit-normalised, but it is a far more convenient convention 
for a supersymmetric theory with superpotential, due to the non-renormalization 
theorem for /’"-terms. Since the superpotential is 


cp 3 

W = Y' 


(4.1) 


it follows that has dimensions °c [mass] 2 / 3 . The Lagrangian must be classically 
scale invariant: in three spacetime dimensions it follows that the Kahler potential 
has dimension 1 and it must scale as K |<3>[ 1 2 . We fix the proportionality constant 
to 

l^! 3/2 i (4-2) 

which gives the kinetic term 

^dn = b K \(t)r 1/2 d(t>d^ 

and the potential 

The Yukawa coupling 

-^Vuk = i0VW + (h.c.) 
does not involve the Kahler potential at all. 

Inclusion of the one-loop renormalization 

Now we apply the RG equations at large |0| at the fixed point, including the small 
one-loop renormalization of the action. (For the running of the Kahler potential, for 


(4.3) 


(4.4) 


(4.5) 


instance, one can use [24], generalizing an earlier formula [25] in the renormalizable 
case.) 

In three dimensions, we have: 

dK 

A — oc A lr 
dA 


where 


2 \W‘ 

m — 


//12 


K 2 , 


2 + m 2 ] , 

(4.6) 

K 2 - 

(4.7) 
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We expand the Kahler potential in powers of \(j> | as 

K(<Z>) = (<P) + (<P) = l( ^\<S>\ i 2 + K W + ■ ■ ■ . (4.8) 

If we assume that the mass scale M = \(j) | 3 / 2 defined by the vev satisfies 


|^| 3 / 2 > A 


(4.9) 


and use the condition <C ~ |0| 3//2 we obtain a self-consistent expansion of 
the logarithm, 


In [A 2 + m 2 ] 


bKA 2 2 | 0|>/2 

, -~-A. , ? ~r 

\(j )\ 3 b K M 


(4.10) 


/, N p- y\2 

We then see ' must scale as kp . Higher terms in the Kahler potential can 
be determined iteratively from the fixed point equation; this generates a series in 
the dimensionless ratio ^A 2 /|0| . All the basic ideas in the determination of the 
A-dependent terms are as in the toy model and the 0(2) model. 


Reduction to Goldstones in the supersymmetric model 

Now we write the leading terms in the large-i effective action in terms of \(j)\ and 
X = arg(0) as in the 0(2 ) model. 


= fc CdflXdj) + v(\<j>\)=b K \<l>\* (dxf + b^— r-+V(|0|) 


I0l : 


101 s 

-I- (higher derivatives) + (fermions) , (4.11) 


with V (101) °c 101 2 . 

For configurations with | (j) | constant, the action is minimized when 

(dx) 2 ~ | (j) | 3 + (higher derivatives) + (fermions). (4.12) 

Eliminating | (j) \ classically and rewriting the effective action in terms of x alone, 
we have a Goldstone effective Lagrangian with exactly the same structure as in the 
case of the 0(2) model, with the exception of the presence of the fermions: 

= b x \dx \ 3 + (lower order in \dx\) + (fermions). (4.13) 


4.3 Decoupling of the fermions 

The fermions, however, decouple from the dynamics. As we show in the Appendix, 
they have a rest energy Eq of order \x\ — 0( x fp), and therefore are heavier than the 
cutoff of our effective theory, and we integrate them out. The fermions obtain their 
large masses from the Yukawa couplings ( |4.5[ ), which allow quanta of R-charge 
carried by the fermions to convert into quanta of R-charge in the % sector. This 
induces a chemical potential for the fermions. 


20 









Importantly, the same Yukawa couplings also impart effective Majorana mass 
terms to the fermions that generate the gap. For most purposes of the large-7 ex¬ 
pansion of the lowest-dimension operators, the key fact is that the low-lying large-7 
sector is described by exactly the same universality class - a conformally-invariant 
effective field theory for the R-Goldstone % - as describes the large-7 sector of the 
critical 0(2) model. 


Rest energies and speeds of the fermions 


The heavy fermions do display a few interesting features at large 7, however. For 
one, their rest mass obeys a precise identity at leading order in 7, 


F n — 2 — - 44 / — -\v\ 

— z dJ R — 2 dJ — 2 IaI j 


(4.14) 


which is necessary for Bose-Fermi degeneracy. The supercharges Q a map from a 
sector of R-charge y to a sector of R-charge y — 1. They do so by removing one 
quantum of R-chargc from the Bose condensate and replacing it with a fermion that 


is almost at rest. Therefore Bose-Fermi degeneracy implies (4.14) to leading order 
in 7. 

The other interesting feature of the fermion dynamics has to do with their prop¬ 
agation speed. Like the bosons, the fermions propagate slower than the speed of 
light. One might have imagined that the fermion and Goldstone speeds are forced 
by supersymmetry to be the same, but this is not so: The fermions actually have 
speed times the speed of light, rather than 4^ as the Goldstones do: 


E f (p)=E 0 ±v f \p\+O(^), 


(4.15) 


where 


Vf = 


(4.16) 


This dispersion relation, particularly the unfamiliar appearance of the negative 
velocity is striking and begs for further explanation. Like the value of the speed 
of the Goldstones, the speed of the fermions is also dictated by superconformal 
symmetry. Consider that the state |7) is the lowest state with R-chargc 37/2 and is 
thus is annihilated by both the energy-raising, R-charge-raising generators Q' and 
the energy-lowering, R-charge-raising superconformal generators S T . However the 
state is non-BPS and therefore can be annihilated by neither S nor Q. Therefore the 
states Q |7) and S |7) are nonvanishing, and have energies A j±\. 

These states are single-fermion excitations on top of a bose condensate of (/)- 
charge 7 — |, so the rest mass Eq = \ \%\ compensates the loss of | of a unit of 

^-charge from the Bose condensate, modulo an error of order J~ 2 . The 0( 1) dif¬ 
ference can thus only come from the kinetic energy of the fermions, which is 


±Vf\p\. On the unit sphere, the role of p in the dispersion relation (4.15) is played 
by the eigenvalue of the Dirac operator, which for the l — \ mode has absolute value 
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1. It follows that Vf can only be equal to \ if superconformal symmetry is respected, 
and that both positive and negative velocities must be present! 

The dynamics of the fermions, particularly the negative velocity modes, is quite 
interesting in its own right. Such modes are widely studied in condensed matter sys¬ 
tems. They are a generic feature of fermions with a chemical potential and Majorana 
mass term, as in our system. The dispersion relation ( A.29| ), for instance, precisely 
appears in famous work by Fu and Kane J26| on the interface between a topological 
insulator and a superconducting material. For a review of such modes and their role 
in condensed matter physics, the reader is referred to p71. 


5 Observables 

5.1 Classical approximation to large-/ states 

We now compute the large-/ expansion of operator dimensions using radial quan¬ 
tization. That is, we use the state-operator correspondence, and the fact that the 
dimension of the lowest operator of charge / is equal to the energy of the lowest 
state of charge / on a unit sphere. 

At leading order, all expectation values of operators are given by evaluation of 
those operators in the lowest-energy classical solution with charge /. By standard 
manipulations in classical mechanics, one may see that the lowest classical solution 
of a system with a global symmetry is always invariant under a "helical" symme¬ 
try - that is, a symmetry under a combined time translation and global symmetry 
transformation. In terms of the Goldstone field, this is simply the statement that the 
lowest-energy classical solution with a given global charge always has the Goldstone 
X varying exactly linearly in time at a rate that is spatially independent. 

This simplifies our search for the classical solution representing the ground state. 
Starting from the Lagrangian density 2z? — b x ( % 2 — ^X 2 ) ■> the conjugate momen¬ 

tum to the field x is 


n 


&X 


= 3b x x(x 2 -(V X ) 2 )' /2 , 


in terms of which the Hamiltonian density reads 


(5.1) 


jfr = 


3V6 


9^(vz) 4 +4n 2 ) 3/2 - 2ih\{yxf +36 n 2 (v x y 


1 1/2 


In terms of X , the Hamiltonian density is given by 


je = b x (r - (v *) 2 ) 1/2 (i% 2 +(v z ) 2 ) 


? _(V X) 2 ^ 4112 + 9*2 ( V*) 4 • 

l ~ 2 6b x 


(5.2) 


(5.3) 


For fixed % or fixed n, the energy is monotonically increasing as a funcion of 
(V^) 2 . It follows that the spherically symmetric state is a global, rather than merely 
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a local, minimum of the energy with fixed total charge. Note that this is a rather 
different situation from that in the system studied in fl7| , where the spherically 
symmetric classical solution for the lowest state at large monopole number is not 
even perturbatively stable. In our case, the spherically symmetric state is not only 
a local minimum but a global minimum of the energy in the configuration space at 
fixed 7. It follows that the lowest state with sufficiently large charge always has 
spin zero in both the 0(2) model and W — <E> 3 model. For a spherically symmetric 
configuration, the Hamiltonian density of the lowest state of the system is 


n = p = 3 b x r 


j 

4 K ’ 


and the classical solution is 




n2 = 


J2 


\fZfb% 4 k 21 Tib x 


(5.4) 


Xo = 



(5.5) 


The total energy of the state is the Hamiltonian density (5.4) integrated over S 


■ 2 . 


E = A (J) 


JT 


^JllKby 


(5.6) 


5.2 Large-/ expansion of the operator dimension 

3 1 

Order J+t and 7 + 2 classical terms 

3 

At order 7 + 2, the energy is just given by the classical energy of the lowest classical 
solution. These contributions come directly from the only two terms in the La- 

grangian with positive 7-scaling, the leading operator @ 1 1 d% 1 3 and Ricci coupling 

_ 2 

(J\ Ric 3 . As shown in Section 3.4 there is no operator of order 7° that contributes 


classically to the energy of the bosonic ground state. So the only order 7° terms 
come from the one-loop vacuum energy, calculated with the Gaussian action for 
fluctuations from the leading-order bosonic terms. We now calculate the value of 
that correction to the vacuum energy. 


Order 7° quantum correction 

We would like to use the leading large-7 terms in the conformal Goldstone action 
to calculate the energy of the lowest state. Since the Goldstone action describes a 
large-7 universality class that is shared by the critical 0(2) model and the W — <t> 3 
model, the following discussion shall apply to both cases together. 

The action for the phase variable x is given by 

■Z’ = b x \dx\ 3 = ^[-^Z^*]+2. (5.7) 
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To quantize it, we break up the field as in section 3.3 i.e.. 


X = Xo + \dXo\ 


(5.8) 


where Xo(t) is the classical solution (5.5) and % is a canonically normalized fluctu¬ 
ation. 

At large-/, we can Taylor expand the action: 


& = b x (dtXdtX - X A 5 2 x) 3/2 

= h\Sx«\ 3 + 3 -^x(-s? + 1a*) i+ o(-§-t) , 


(5.9) 


where we have dropped a total derivative term proportional to d z %. 

Wick rotating t —> iz and computing the determinant, we find that the one-loop 
correction to the vacuum energy is given by the usual Coleman-Weinberg formula 
applied to a minimally coupled massless boson with propagation speed times the 
speed of light: 


2 ^ logdet f-<9 T 2 -^A 5 2 ) = 


i oo 

—= £(2m)\W+i). (5.io) 

/v / £=0 


T /2 

(We have let T = f _J r p dx denote the total extent of (Euclidean) time.) By using 
the ^-function regularization, the sum can be renormalized to be 


£(2f+l)V^+l)= £(2m)vW+l) 

£=0 £=\ 

1 
4 


= £ l2f + 2£+j + 0(r 2 ) 

i=\ v 

2£(—2) + 2£(—1) + -£(0)+ 


°° / 1 
+ £ (2£+l)^/£(£+l) -2£ 2 -2£- - 


(5.11) 


l=\ 

= 0.01508 = -0.30675. 
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Multiplying by A__, we h ave the one-loop renormalized vacuum energy of the 
bosonic 0(2 ) model in the large-/ sector: 


A(/) 


= -0.108451 . 


(5.12) 


O(J 0 ) term 

Thus the energy of the lowest state has the asymptotic expansion 
A(/) = c, 3/ + 2 +c,i/ + 2 -0.108451 + 0(/-2) . 


(5.13) 
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It is sometimes helpful to re-express the structure of the asymptotic expansion 
in terms of a sum rule. We therefore characterize the expansion for the ground state 
energy (|5.13[) as the sum rule: 


j2a ( j )-(j + J 4 + Js) + A (/+ l ) = 0.04067 + 0 ( 7 -!). 

( 5 . 14 ) 

As we have emphasized throughout, there are no operators scaling as 7°, so the order 
7° term is universal. 


Noncontribution of higher-loop diagrams at order 7 


The Casimir contribution to the energy is the entire contribution at this order. It is 
easy to see that there are no higher-order diagrams contributing to the energy at order 
7°. Once we write the action in terms of the canonically normalized fluctuation field 

3 3 

X, we see that every cubic vertex in a Feynman diagram scales as \r)Xo\ 5 p * 
at the largest, and every quartic vertex as \dxo[ 3 p~ 2 at the largest. Therefore, a 

_ 3 

two-loop vacuum diagram contributes as 7 2 and smaller. 

Corrections to the ground state energy from subleading explicit terms ( |3.36[ ) are 
larger but still too small to contribute at order 7°: The operator &\ = Ric3|7^ gives 


(Ric) 


(dX ) 2 


oc p 1 when expanded to quadratic order in fluctuations. It affects the 


dispersion relation and Casimir energy at order 7 1 . 


5.3 Energies of excited states 


It is now simple to work out the energy of the excited states and their conformal 
representation theory, i.e., which states are primaries and which are descendants. 


From ( |5.9[ ), we see that at leading order the Lagrangian is just that of the free fields, 
so the equation of motion is just 


1 2 

X,tt = + 2 V " Z ' 


(5.15) 


Note that the dispersion relation is nonrelativistic; the propagation speed of the 
Goldstone boson is times the speed of light. We will see momentarily that this 
value is the only one consistent with the spontaneously broken conformal symmetry. 
The dispersion relation in terms of frequencies and spins is as usual on a sphere: 


co e =^£(£+l). (5.16) 

The energies of the excited states are therefore the energy of the ground state, plus 
the sum of a set of frequencies (Og. 

We start by examining the first-excited states. Take l — 1 for example. The first 
excited states have one excitation of spin 1 that increases the operator dimension by 
1, which corresponds to acting with d on the primary operator to make a descendant. 
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The speed of the Goldstone, 1 / y/2 plays an essential role here because otherwise 
there would be no spin-1 oscillator which would increase the anomalous dimension 
by 1. The conformal raising operator P' fl , in other words, is just the Goldstone mode 
at l — 1. Therefore the criterion for a state to be a conformal primary is that it has 
no £ — 1 modes excited, with only the modes of £ — 2 or greater occupied. 

We can now compute the weights of primaries. Acting with a single £ — 2 os¬ 
cillator gives an energy of A(/) + a/ 3, which is the lowest primary state above the 
large-/ vacuum. This state must be primary, because there is no state of charge / 
with energy A(/) + \/3 — 1. We can construct scalar primary states as well. The 
lowest-energy way to do this is by taking the singlet in the symmetric product of 
two £ — 2 representations. Therefore the lowest excited scalar primary is obtained 
by acting with two £ — 2 oscillators, and therefore has energy A(/) + 2 a/3. The 
lowest vector primary is obtained by acting with one £ — 3 oscillator and one £ — 2 
oscillator, and has energy A(/) + y/3 + y/6. 

All these states have subheading large-/ corrections coming from loop effects 
and explicit operator corrections in the Lagrangian, of course. These contribute to 
E — A (/) at order / 2 and / 1 , respectively. 


6 Large-J analysis of some other models 

6.1 4D theories and their anomaly terms 

Consider a CFT in four dimensions with a U( \ ) global symmetry, in which Weyl 
anomalies are present in the underlying CFT. And suppose further that the large / 
universality class of this CFT is described again by the a four-dimensional version 
of the conformal Goldstone system studied in the previous sections^ In such a 
CFT, the anomalies must express themselves consistently at the level of the large-/ 
effective dynamics. The result is that the anomaly coefficients control certain terms 
in the large-/ expansion of the CFT. In four dimensions, the a-anomaly directly 
dictates the coefficients of certain terms in the large-/ effective action in flat space; 
certain curvature-dependent terms of order /° in the large-/ effective action are also 
dictated by the Weyl anomaly coefficients, of both a- and c-type. 

The analysis of the Wess-Zumino term for anomalous Weyl symmetry has been 


12 Though we do not know of any controlled example of such a CFT, we can certainly imagine they 
may exist. For instance, an asymptotically free SU (2) gauge theory with two pairs of complex 
scalar fields, each transforming as a doublet, may have a phase transition as the bare mass-squared 
parameter is varied from positive to negative and the U( 1) global symmetry acting on the doublets 
goes from unbroken to spontaneously broken. Assuming such a phase structure exists, then the 
boundary between the two presumably has a description in terms of a four-dimensional version of 
the conformal Goldstone theories we have discussed in previous sections. 
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carried out in [28] (see their Eq. (2.8)) and the result is that 


S W z[ggv,X',a,c] 



— a 


ln(\dx\)E 4 (g) 


+ 4(R^(g) - l2g pv R(g))^~ 2 8^D.8 v D. 
+ 4fr 3 (d9) 2 Da-2fr 4 (d9) 4 


+ c ln(|<9^|) l¥ 2 (g) j, 


( 6 . 1 ) 


where g is the background metric, £4 and W^ v p G are respectively the Euler density 
and Weyl tensor, and for us, 

n=\dx\. ( 6 . 2 ) 

The point here is that the logarithm of \d%\ plays exactly the same role as the 
effective dilaton T in the Wess-Zumino anomaly term, that the logarithm of the 
overall scaling modulus plays in the same formula when there is a supersymmet¬ 
ric moduli space of vacua. It would be interesting to learn from this substitution 
whether any interesting insights can be gained from unitarity constraints on Gold- 
stone scattering, perhaps generalizing the recent celebrated proof of the a-theorem 
in four dimensions [29]. 


6.2 Monopoles and S-duality in three dimensions 

Dualities at large J Like many known theories of a compact boson in three dimen¬ 
sions, the conformal Goldstone system has a duality transformation to an Abelian 
gauge theory in three dimensions, where the Noether current maps to the monopole 
current, i.e., the Hodge dual of the field strength tensor, J? p ^ E pva [d A) va . Then 
the total Noether charge of a state maps to the magnetic flux on the sphere. In the 
language of operators, the Noether charge of an operator in the original variables, is 
the monopole number of the operator. 

We can use coordinate- and Weyl-transformation properties in order to see how 
the field strength tensor can be expressed in terms of X■ We take here a convention 
where £012 = 1 and thus has Weyl weight 0. Also as F^ v has Weyl weight 0 and 
we are in three dimensions, so that \f\g\ has weight —3. At leading order in the 
derivative expansion, we find that Weyl-invariance, diffeomorphism covariance, and 
charge quantization uniquely determine the relation: 

Fgv — V2\d X \ (*dx)gv — \^x\ Xi (6.3) 

and the inverse relation is: 

3nx = 2j ifT l/2 (*n. = ^ i*r 1/2 \/i*lw-f vp > < 6 - 4 ) 

where the numerical factors have been chosen so that 

\F\ 2 = \dx\ 4 . (6.5) 
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The duality means that the effective Lagrangian for the field strength is immedi¬ 
ately derived from the leading Goldstone action in Eq.(|3.29[): 


^ = b y \F | 3/2 + .... 


( 6 . 6 ) 


This is consistent with the fact that the Weyl weight of the Lagrangian is 3. An 
immediate consequence of the unusual form of the action is that the dimension of 
the lowest-lying monopole operator, scales as monopole number to the |, for large 
monopole charge p] 

Note that the relation between and F^ v is exact, but the relationship between 
and % depends on higher terms in the Lagrangian. However the modifica¬ 
tions have a sub-leading effect at large magnetic flux number: corrections are sup¬ 
pressed by inverse powers of constant field strength, with numerators proportional 
to derivatives of the field strength. 


7 Conclusions 

We have performed a renormalization group analysis proving that certain simple 
bosonic and supersymmetric systems are described at large charge density by a sim¬ 
ple conformal Lagrangian for Goldstone fields. In the limit where the charge density 
is large compared to the infrared energy scale, the system is weakly coupled both 
classically and quantum mechanically, and can be quantized straightforwardly in 
perturbation theory. In the limit of large charge 7, the leading large-7 expressions 
for all quantities, such as energies and anomalous dimensions, are controlled by 
leading terms in the effective action. We do not know at present how to calculate 
the coefficients of those terms analytically from first principles, and we expect them 
to differ between conformal field theories. Lor instance, the value of b describing 
the large-charge sector of the W — <f> 3 supersymmetric model likely differs from the 
value of b describing the large-charge sector of the three-dimensional XY model. In 
even dimensions certain coefficients in the large-7 effective action can be computed 
in terms of some intrinsic data of the CFT such as anomaly coefficients. 

We have also computed the dimensions of excited primary states, and found that 
they are given up to order 7° by energies of free oscillators in the A > 2 spherical 
harmonics on S 2 . 

There are a number of interesting questions to investigate in the future. 

Since some features of the three-dimensional model may appear counterintuitive, 
it would be interesting to compare them with the analogous properties of various 
known conformal models in two dimensions. The complete solvability of these 
models would give more confidence in the consistency of our framework. 

We may hope our framework is powerful enough to provide insights in the large- 
7 behavior of other strongly coupled CFTs which are in general not tractable with 

13 The same scaling has apparently been derived for operators of large monopole number in three- 
dimensional fixed points flowing from weakly-coupled gauge theory or Chern-Simons matter theory. 
We thank Ethan Dyer for communicating this result to us. |2TJ . 
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known methods. Natural examples are 2D cr-models with degenerating cycles, non- 
supersymmetric Chern-Simons-matter theory at finite rank and level [18] 191, and 
the (2,0) theory at finite N. 

A brief list of further interesting directions is given below. 


Constraints on large-2 Lagrangian parameters. In some cases, one can con¬ 
strain terms in the large-7 effective theory in terms of the microscopic theory. This 
is principally the case for four-dimensional theories, where certain terms are dic¬ 
tated by anomaly coefficients. It would be valuable to understand whether this can 
be done more generally by matching correlation functions. It would be nice, e.g. 
to match the |x|~ 4 term of a current two-point function in three dimensions with a 
current correlation function in the large-7 effective theory, in order to calculate the 
/^-coefficient in terms of CFT data. 


Theories with moduli spaces and bps operators. For theories with moduli spaces 
and (equivalently) an infinite ring of BPS primary scalar operators, our methods cer¬ 
tainly cannot add anything to the prediction of the leading-order energies of the 
BPS states, since their dimensions are dictated by an exact formula to be equal to 
the //-charge, A = \Jr\. However large-7 methods can certainly shed light on re¬ 
lated questions not controlled directly by superalgebraic considerations, such as the 
anomalous dimensions of states lying just above the BPS bouncf* 4 ] 


Comparison with other results on the critical 0(2) model. One would like to 
see how our asymptotic formulae for the operator dimension may fit with other 
approaches to understanding the 0(2) model at the critical point. For one, the model 
can be simulated on a lattice and our results tested numerically. Also, much recent 
progress has been made on the conformal bootstrap both for the three-dimensional 
critical 0(N) models and the 0(2) model, specifically [ 131, as well as for the W — 
<f> 3 theory [3lJ 32|, and it would be desirable to see how such analytic bootstrap 
methods may confirm and/or complement the results obtained in the present paper. 
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A Fermions and supersymmetry at large R-charge 


In this section of the Appendix, we work out some general features of the situation 
in which we have a theory with four supercharges and an exact R-symmetry, and 
we examine the lowest state at large R-charge and the energies of the fermionic 
excitations above it. The main focus of this section is to derive at the Lagrangian 
level the large gap in the fermion sector in the sector of large R-charge. 


A.l Quantum mechanics with four supercharges 


We study a quantum mechanical theory with two complex supercharges Qa,Qce , 
where we take a to run over the indices a e {t-!}. The SUSY algebra is 

{Qa,Qp} = 0 {Ql,Ql} = Q {Q a ,Q\} = 2H8 a p. (A.l) 

For a chiral multiplet ((f). (//. F), the SUSY transformations are 


Qa'fy — 

Gi-0 = 0 

Qa-Yp = V2e ap F 

Qa ' = 0 > 

QaF = 0 

Qi -F = -iV2 e a p y/p 


Qa ' 0 — 0 j 

Qlc-<l> = -iy /2 
Qay/'l = V28 a p$ , 

Qa~¥p = ^£apF 
Qa-F = ~iV2£ a p ¥p 
Ql-F = 0 . 


The simplest model with this super algebra and field content is the 0+1-dimensional 
Wess-Zumino model, which is defined by a Kahler potential K((j),(j)) and holomor- 
phic superpotential W. Taking the Kahler potential to be the flat, unit-normalized 
one (j>) = we have 


Fk = -4 2; • Q\ ■ Q\ ■ K{ (j) , (j) ), = | <j) | 2 + |F \ 2 + i 1/4 Ya + (total derivative), 


(A.3) 

L w = ^Q r Q i -W^) = -W\8)F + i W"(0) (A.4) 

% = l^ = -^Gj-Gj-W(0)=-W / (0)F+iW // (0)v^v|- (A.5) 

The total Lagrangian is L = L K + L W + L w = Lb os + Lferm, where 

T bos = |0| 2 + \F\ 2 -FW\8)-FW\§) (A.6) 

£ferm = l'v4^a + 1 ^"(0) 1^; + iW" ($) V'j• (A.7) 

After integrating out the auxiliary fields, we have 

F = W , (0), F = W'{<$>), (A.8) 
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L hos = \<j)\ 2 -\W\<l))\ 2 . 

Now specialize to the case 




gl -1 

(q+l)^q 


ty q+X ■ 


(A.9) 

(A. 10) 


where we have chosen our normalization for g to simplify the formulae. For the 
superpotential (]A. 10[) 


F 


g c i- ] ^ 

Vq 



(A. 11) 


^bos = |0| 2 -^|0| 2? (A. 12) 

<T 

^ferm = \j/+ i ^/q g q ~ 1 ^ VW + ^VW 1 f* -1 vjtyf- (A.13) 

Writing our complex scalar (j> in polar coordinates as in (|3.3[), we have 


p 2q-2 

Lbos=a 2 X 2 - a 2q 


q 


^ferm — / \j/ ' \j/ i y ^q g^ Cl 


9-ln9-l 


(A. 14) 


exp{/ (q - l)x} YtVi + exp{-i (q - 1) x} vf w\ 

(A. 15) 


We now make an additional redefinition to eliminate the nonderivative coupling of 
the R-Goldstone x to the fermions: 

V = exp{-^^ x} ¥ , ^expj+^A^zlV' 4 , 

(A. 16) 

y/ = ex p{- l -^ 2 ^x}V, t// 4 = exp{+^=^z}Y rt , 

in terms of which the fermonic Lagrangian is 

Lf ena = iY^+^^XV^Y + iVqg q ~ laq ~ 1 (Vm + V^^-l) ' (A ' 17) 

For fixed X, the frequency of the a-oscillator scales as X- so we would like to in¬ 
tegrate it out it of the system and obtain an effective theory of the other degrees of 
freedom. Working classically, we integrate out a by eliminating it classically from 
the Lagrangian by minimizing its energy at fixed x, giving 

(ga) 2q ~ 2 = X 2 , a = g~ 1 \%\& (A-18) 

so for states with a in its ground state we have 

1 <7 — 1 2g 

L b os = ^~ -Izb- 1 , (A. 19) 

g 2 q 

L fsIm = i¥ t P+^ ] -X¥ t ¥ + iVq\X\ + ’ ( A -20) 
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which we now write as 


Lferm = l ^ + JU \j/ J + lAf 


t, 


W]Wi + w\wl 


(A.21) 


Af = v^lATl » (A. 22) 

where /i is a kind of pseudo-chemical potential for the y/-chargc and y/ is a sym¬ 
metry breaking mass term that violates t/r-charge. The parameter fl can only be 
understood as a true chemical potential in the limit M —> 0. For fixed q we are not 
within the confines of such a limit, so /i does not admit a consistent interpretation 
as a chemical potential in this context. 

The energy of the fermionic excitations is 

£= VV 2 +M 2 = 2+i \%\ . (A.23) 

There are two degenerate modes with this energy, one with each possible eigenvalue 
±1 of the internal 517 (2) magnetic quantum number. Notice that the excitation 
energy of the fermion is parametrically large as a function of \x\. So if the upper 
limit on our energy is some fixed scale A independent of J, then we are justified in 
treating the fermionic degrees of freedom as frozen in their ground states. 


A.2 Lift to field theory 


Now we lift the model to 3 + 1 dimensions; this is just the four-dimensional Wess- 
Zumino model of a single chiral superfield a homogeneous superpotential of the 


form (A. 10). This model has a Landau pole for q = 2 and is non-renormalizable for 
q> 3, but at present we are only considering semiclassical aspects of the theory, and 
we can define the model quantum mechanically with a cutoff A if desired. 

The form of the Lagrangian density is exactly the same as that of the Lagrangian 


in the previous Section A. 1, except that the fermion kinetic term acquires a piece 
with spatial derivatives 


JZbos = 

^fern 


\wm 


vwi + iw"^) v' t Vf 


(A.24) 
(A.25) 


Specifying to the superpotential (A.10), going to polar field variables as in (3.3), 


and defining hatted fermions as in (A. 16), we have 


Lferm = ~i d, L \jf - W + iM 


/V y\ . A T A T 

W]Wi + W\wl 


(A.26) 


q -1 




M = \fqg q 


V- 


(A.27) 
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Now restrict to the lowest-energy state with given /^-charge, for which the gradi¬ 
ent of x is constant and purely in the time direction, and a is constant as well, satisfy¬ 
ing (A. 18). Then letting p e (0,A) and using xfr = i//’T° and P 4 F° = —r°r A = + o A 
when acting on Weyl spinors of a definite chirality, we have 


^ferrn = ify'+ i^C A d A fy + jU ij/Ujf+iM ( ijf^ + fi 


(A.28) 


with M and p given as before by (A.22). The excitation energies of the fermions are 
then 


E± = y/(\ P \±p) 2 +M*, 


\p\ = yJPAPA ■ 


(A.29) 


The energy of a fermion at rest in the reference frame of the charge density, has 
energy 


So =^1*1 • 

For small momenta \p\ <C \%\, the dispersion relation is 


E± ~ Eo±v/\p\ , 


v f 


q- 1 

< 7+1 


(A.30) 


(A.31) 


The same dispersion relation ( A.29[ ) holds when the theory is reduced to any 
lower dimension D, with spatial momenta taking values in D — 1 dimensional vec¬ 


tors. On reduction to 0 + 1 dimensions we recover (A.23). 


A.3 Nontrivial Kahler potential 

It is clear that the parameters M, /i in the dispersion relation, expressed in terms of %, 
are independent of the coupling g. It follows immediately that they are independent 
of the normalization of the flat Kahler metric as well, because the normalization of 
the flat Kahler metric can always be absorbed into a redefinition of g, via a rescaling 
of the fields of the multiplet (0, t j/,F). 


It is slightly less obvious, but still easy to see, that the value ( |A.30[ ) still holds 
exactly, even if the Kahler potential has a general homogeneous form K °c \<t>\ 2a - By 
redefining 


0 = 0 ' 


(A.32) 


-q +1 


we obtain a quadratic Kahler potential K °c \(j)\ 2 and a superpotential W 0+“. 
We then recover the same system as in sections ( |A.1[ ),( A22| ), with % replaced by the 
phase % of (j), and with q replaced with 


<7ne w = 2 ± i -l- 


(A.33) 


All formulae for the fermion dynamics from section A.2 then hold, with the value 


of a replaced by q new . In the case of interest to us in the present paper, namely the 
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three-dimensional superconformal fixed point, the exponent in the Kahler p otentia l 
is a — |, and the superpotential is cubic, so in our notation q — 2. Thus using (A.33), 
we have <7 new = 3. 

We therefore have 


A i = Z, M = (A.34) 

Eq — 2\x\ , v f = 1 -. (A.35) 


A.4 A general supersymmetric identity on the rest mass 

The formula for the energy Eq in ( |A.30[ ) is directly dictated by supersymmetry. Even 
though all supercharges are realized nonlinearly, there are no massless goldstini in 
the system. This is because the supercharges do not act on the fixed-i Hilbert space; 
rather, they map from the Hilbert space with one value of the /^-charge to another. As 
a result, the rest mass of the fermion is given not by zero, but by the change in energy 
under the removal of a unit of 0-charge. Since the ratio between the normalisation 
of 0-charge to that of /^-charge is ^y, the condition of Bose-Fermi degeneracy is 


dA(J) _ q+\dA(J) 
cIJr 2 dJ 


(A.36) 


so that the energy of a zero-momentum fermion can exactly compensate the loss of 
one unit of ^-charge from the Bose condensate. 
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